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ABSTRACT 

Context. We explore the dynamical evolution of a planet embedded in a disk surrounding a star part of a binary system where the 
orbital plane of the binary is significantly tilted respect to the initial disk plane. 

Aims. Our aim is to test whether the planet remains within the disk and continues to migrate towards the star in a Type l/Tl mode in 
spite of the secular perturbations of the companion star. This would explain observed exoplanets with significant inclination respect 
to the equatorial plane of their host star. 

Methods. We have used two different SPH codes, vine and phantom, to model the evolution of a system star+disk+planet and com¬ 
panion star with time. 

Results. After an initial coupled evolution, the inclination of the disk and that of the planet begin to differ significantly. The period 
of oscillation of the disk inclination, respect to the initial plane, is shorter than that of the planet which evolves independently after 
about lO'^ yr following a perturbed N-body behavior. However, the planet keeps migrating towards the star because during its orbital 
motion it crosses the disk plane and the friction with the gas causes angular momentum loss. 

Conclusions. Disk and planet in a significantly inclined binary system are not dynamically coupled for small binary separations but 
evolve almost independently. The planet abandons the disk and, due to the onset of a significant mutual inclination, it interacts with 
the gas only when its orbit intersects the disk plane. The drift of the planet towards the star is not due to type l/Il with the planet 
embedded in the disk but to the friction with the gas during the disk crossing. 

Key words. Protoplanetary disks — Methods: numerical — Planets and satellites: formation 


1. Introduction 

According to observational surveys, about half of solar-type 
stars reside in multiple stellar systems ( |Raghavan et al.|[2010| l 
with the frequency declining to roughly about 30% for less mas¬ 
sive M stars ( |Lada 2006| l. This frequency is suspected to be 
higher among young stars ( |Reipurth]|2000[ ) and the subsequent 
decay to the present fraction may be due to dynamical instabil¬ 
ity or gravitational encounters with other stars. The presence of 
circumstellar disks around both components of binaries does not 
seem to be strongly affected by the gravitational perturbations of 
the companion. Spatially resolved observations of disks in bina¬ 
ries in the Orion nebula cluster ( [Daemgen et al. |2012| l suggest 
that the fraction of circumstellar disks around individual compo¬ 
nents of binary systems is about 40%, only slightly lower than 
that for single stars (roughly 50%). On one side the secondary 
star perturbations may lead to disk truncation (|Artymowicz & 


Lubow 1994| 1996||Bate 2000|l reducing its lifespan, on the other 


side the presence of a circumbinary disk may fed the truncated 
circumstellar disks through gas streams so that their dissipation 
timescale does not vary significantly compared to single stars 
( |Monin et al.|2007) and they have a chance to form planets. 

Up to date, 97 exo planetary systems have been discovered 


around multiple stars (Rein 


.ys i 

2012 n triggering statistical com¬ 


parisons with planets around single stars (|Desidera & Barbieri| 


* github. com/OpenExoplanetCatalogue/open_exoplanet_ 
catalogue/ 


|2007[ [Roell et^|2012[ ). The analysis indicates that planetary 
masses increase for smaller stellar separation while for wide bi¬ 
naries the physical and orbital parameters appear very similar. In 
effect, planet formation in close binary configurations is a com¬ 
plex problem (see Thebault & Haghighipour|2014| for a review) 
while it is expected that for large separations the effects of the 
binary companion on the planet growth may be less relevant. 
However, if the orbit of the companion star is significantly in¬ 
clined respect to the initial protostellar disk plane, which is pos¬ 
sibly co-planar with the star equator, this may lead to important 
dynamical consequences on the final orbit of the planet, in partic¬ 
ular its inclination respect to the star equator. It has in fact been 
suggested ( |Batygin||20I2| l that the evolution of the protoplane¬ 
tary disk under the perturbations of the binary may be responsi¬ 
ble for the observed spin-orbit misalignment of some exoplan¬ 
ets. According to |Triaud et al.| ( [2010[ ) and [Albrecht et al.] ( |2012| l 
about 40% of hot Jupiters have orbits significantly tilted respect 
to the equatorial plane of the star. A fundamental requirement of 
the model of |B aty gin] ( |20 1 2| l is that the planet resides within the 
disk during its evolution. In this way it would continue its mi¬ 
gration by tidal interaction with the disk and, at the same time, 
it will follow the inclination evolution of the disk. Once the disk 
is dissipated its migration will stop, while the evolution of the 
inclination will continue, though in a different fashion, until the 
binary companion is possibly stripped away leaving the planet 
on an inclined orbit. 
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Previous studies seem to suggest that indeed a giant planet is 
forced to evolve within its birth disk even in presence of exter¬ 
nal perturbing forces. Any relative inclination between the planet 
and the disk plane, induced either by a planet-planet scattering 
event or by resonances (Thommes & Lissauer|2003|, is quickly 


damped by the disk leading to a realignment of the planet orbit 


with the disk plane (jCresswell et al.| 

20071 Marzari & Nelson] 

20091 Bitsch & Kley|2011[ Rein|2012 

[Teyssandier et al.|2013|l. 


However, the precession rate of the planetary orbit due to the in¬ 
teraction with the binary star is faster than that of the disk. As 
a consequence, if the damping of the disk on the planet orbit is 
not strong enough to keep it within the disk plane, the planet will 
evolve independently under the secular perturbations of the sec¬ 
ondary star and develop a significant relative inclination respect 
to the disk. The assumed type II migration will not occur in this 
case since the planet is not embedded in the disk and cannot open 

a gap. _ _ 

In a recent paper, Xiang-Gruess & Papaloizou (|20I4|l have 


investigated the evolution of a disk and massive planet under 
the influence of an inclined binary companion. According to 
their SPH simulations, the planet and disk maintain approximate 
coplanarity during the evolution of the system and the planet 
would quickly migrate towards the star and stop on a close orbit. 


This findings seem to confirm the model described in Batygin 
( [2012| l for small separations between the binary components. To 
further test these findings, we have performed numerical simu¬ 
lations with a setup similar to that adopted in |Xiang-Gruess & 


Papaloizou 

mu 

1 and used two different SPH codes: vine (Wet- 

zstein et al. 

2000 

Nelson et al.|2009b and phantom ([Price & Fed- 

errath|2010 

Lodato & Price 201 Ob. Our runs extend over a much 


longer timescale and we adopted a higher resolution. The simu¬ 
lations show that the initial coplanarity is maintained only for a 
limited amount of time and that the planet definitively detaches 
from the disk plane as soon as the disk completes a quarter of its 
precession period. The secular perturbations of the binary com¬ 
panion in an inclined orbit overcome the damping force of the 
disk and the planet evolves independently of the disk. A signifi¬ 
cant mutual inclination develops between the disk and the planet 
orbital plane but this does not halt the orbital migration. When 
the planet crosses the disk it is still affected by the gas via dy¬ 
namical friction and its inward drift continues even if at a slower 
rate compared to that computed by |Xiang-Gruess & Papaloizou| 
( |2014| l in the initial stages of the system evolution. 

In Sec. 1^ we study the theoretical model of the precession 
rates in a star-i-disk-i-planet-star system and compare evolution 
timescales. In Sec. [^we describe in detail the model setup of 
the numerical simulations. Then, in Secs. andwe show the 
results of the long term high resolution models for different bi¬ 
nary inclinations, and in Sec.j^we discuss our results and their 
implications. 


2. Relevant timescales 

We compare in this section the timescales of precession of the 
angular momentum of the planet computed within a pure N- 
body problem and that estimated for the disk. This may give 
important hints on the forces that may try to separate the planet 
from the disk and give a better insight on the results of our SPH 
simulations. In computing the precession timescales, we con¬ 
sider a planet of mass Mp initially on a circular orbit around a 
star of mass M* at a semimajor axis a, with a distant binary com¬ 
panion of mass Mb, semimajor axis a\,, inclination /b, and eccen¬ 
tricity Cb which we assume to be 0 at the beginning. The pertur¬ 
bations of the companion star are strong in particular for incli- 
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nations higher than /b ~ 39° when the Kozai-Lidov mechanism 
( |Kozai|[r962[ |Lidov|[T962[ ) forces wide oscillations of both the 
inclination and eccentricity of the planet whose values are tied 
by the conservation of - ^(1 - e^)cos(/p). The maximum 

eccentricity reachable by the planet is e^ax - -^1 - 5/3 cos^ . 
However, short range forces like General Relativity and tidal dis¬ 
tortions (which we neglect here), and the dumping effect of the 
protoplanetary disk, tend to reduce this value. 

The Kozai cycles occur at a characteristic period given by 
( [Kiseleva et al.|1998] ) 


Pi 

Pkz - 

P 




M*+Mb 


Mb 


QpU/ \Mb/' 




( 1 ) 


(see also Ford et al. 20001 where Pp and are the orbital periods 
of the planet and binary. The dynamical evolution of the planet 
can also be viewed from the point of view of the orbital angular 
momentum. According to Storch et al. ( |2014| l, the planetary or¬ 
bital angular momentum vector precesses around the binary axis 
(Lb) at a rate that, in the absence of tidal dissipation, is approxi¬ 
mately given by 
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( 2 ) 


where eo is the initial eccentricity of the planet, ffp^ is the incli¬ 
nation of the planetary orbit respect to the binary orbit, and 
is its initial value. The precession of the orbital angular momen¬ 
tum approximately translates into an oscillation of the planet in¬ 
clination respect to the initial plane which is the disk plane. The 
timescales estimated by eq.[T]and|^are then comparable. 

Also the disk is affected by the binary companion gravity and 
the disk axis Ld precesses around the binary axis Lb, assuming 
that it satisfies the condition for a coherent behavior as a solid 
body, with a period given by 



1 

K cos 6>db 


^ Itt 


1/2 





1 

3/8cos0db 


(3) 


(see jBate et al.|2000l|Lai|20r4l l, where rout and Rout are the radius 
and Keplerian period of the outer disk edge respectively, and 0db 
is the inclination of the disk respect to the orbital plane of the 
binary. The relative precession rate is given by 



CM* 




cos ^db 


(4) 


Even in this case, the precession can be translated in a peri¬ 
odic oscillation of the disk inclination respect to the initial plane 
which can be assumed to have been the equatorial plane of the 
primary star. 

If we now compare the two precession rates, that of the 
planet and that of the disk, we And that the initial ratio between 
them, under the assumption that the planet formed in the mid¬ 
plane of the disk (0pb = - 0db) in a circular orbit {e - - 0), 
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is given by 


pprec 

_ 

ppi-ec 


fipb 


^AnP-^ 
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(5) 


where in the last step we assumed that the binary truncates by 
tidal interaction the disk maximum radius to 1 /3 of the binary 
separation, which is a good approximation for binary separations 
less than ~ 300 au ( Artymowicz & Lubow| 1994 1. 

For t + to and eo = 0 the ratio is evolving as 


pT cosflf; 

d 2- 


I. 


cos0db\ab/3 




1 ^ sin^ 61“ 


I -e^ I sin^flib 


(6) 


If we consider a typical case of an equal mass binary made 
of solar-type stars on a circular orbit and separated by 100 au 
and a planet orbiting around the primary with a = 5 au and em¬ 
bedded in a disk with outer radius of 30 au, from eq. |^the two 
period differs by about a factor ten. This is a large difference and 
the disk damping must overcome the tendency of the planet to 
evolve on a slower timescale and then out of the disk plane. 


3. Model description 

In this section we briefly describe the two numerical algorithms 
we have used to model the system planet-i-disk-i-binary compan¬ 
ion and outline their differences which might affect long term 
results. 


3.1. SPH codes and model setup 

To model the evolution of a Jupiter-size planet embedded in the 
disk surrounding the primary star of an inclined binary system 
we have used two different SPH codes: 

- VINE ( |Wetzstein et al.|[2009[ [Nelson et al.||2009l l, which is a 

hybrid N-body/SPH code, updated to improve momentum 
and energy conservation as described in [Picogna & Marzari 
( |20T3] |, _ 

- PHANTOM ( [Price & Federrath|2010[ [Lodato & Price|2010| l, a 
modem SPH code which models the massive bodies as sink 
particles ( [Bate et al.|1995| ). 

Both these codes solve the hydrodynamic equations in the La- 
grangian approach by replacing the fluid with a set of particles 
(see [Price|2012[ for a review). 


Equation of state — A locally isothermal equation of state, sim¬ 
ilar to that described in [Peplihski et al.[ ( [2008| l, is adopted in all 
simulations 


Cs = 


hsrshpfp 


[(h,rfi" + (/ipTp)"]!/" 




'2 H- 


(7) 


where Tj, rp are the distances of the fluid element from the pri¬ 
mary star and the planet, respectively, and h^, hp are the circum- 
stellar and circumplanetary disk aspect ratios, and f2p are the 
circular Keplerian orbits of a fluid element around the star and 
planet, and n - 3.5 is a non-dimensional parameter chosen to 
smoothly join the equation of state near the planet (cs = hprpQ.p 
for Tp << rf) with that of a flat circumstellar disk (cs = h^rsflf). 


For the simulations in this work we choose a disk aspect ratio 
of hs = 0.037, and a circumplanetary scale height of /zs = 0.6. 


Accretion — vine treats the stars and planet as N-bodies. The 
mass accretion has been modelled for the stars inside a radius 
of 0.5 au, but not for the planet. 

In the simulations performed with the phantom code, the 
planet (and stars) are all modelled as Lagrangian sink particles 
( [Bate et ak 1995| l. A sink particle is evolved as an SPH parti¬ 
cle, but it experiences only the gravitational force. A gas particle 
that comes within its accretion radius, - 0.075 au, can be 
accreted if 


- it is inside the Hill sphere of the sink particle, 

- its specific angular momentum is less than that required to 
form a circular orbit at the accretion radius, 

- it is more bound to the candidate sink particle than to any 
other sink particle. 


Viscosity — An artificial viscosity term is introduced in SPH 
codes in order to correctly model shock waves that inject en¬ 
tropy into the flow over distances that are much shorter than a 
smoothing length, and to simulate the evolution of viscous disks. 
The term broadens the shock over a small number of smoothing 
lengths and correctly resolves it ensuring at the same time that 
the Rankine-Hugoniot equations are satisfied. In this way it pre¬ 
vents discontinuities in entropy, pressure, density, and velocity 
fields. 

The implementation of the viscosity term is slightly different 
in the two numerical codes, vine adopts the standard formula¬ 
tion introduced by [Monaghan & Gingold[ ( [1983P where, for ap¬ 
proaching particles (Va* ■ Ta* < 0), an artificial viscosity term is 
introduced in the momentum and energy equations 


Hfli - 


(-ac„i,iAai,+Pnlp 

Pab 

0 


Vafo ■ fafc < 0 

Vfl* • faZ, > 0 


( 8 ) 


where pat is a velocity divergence term 


Pab 


byab ■ tCab 
+ 


(9) 


ab 


with - 0.01. This implementation conserves total linear and 
angular momentum and vanishes for rigid body rotation, a is the 
linear term (bulk viscosity) and it dissipates kinetic energy as 
particles approach each other to reduce subsonic velocity oscilla¬ 
tions following a shock, fi is the quadratic term (von Neumann- 
Richtmyer like viscosity), which converts kinetic energy to ther¬ 
mal energy preventing particle mutual penetration in shocks. 

PHANTOM uses a more general formulation of dissipative 
terms (Monaghan |1997| l. It is built on an analogy with Rie- 
mann solvers, where the dissipative terms of conservative vari¬ 
ables (density, specific momentum, and energy), that experience 
a jump across a shock front are multiplied by eigenvalues similar 
to signal velocities (fsig)- The viscosity term for the momentum 
equation becomes 


1 CrV^igyah ■ ^ab 

Hfli = -- 

2 Pab 

where the signal velocity is given by 


( 10 ) 


_ jCs,a 3" Cs^b fyah ' ^ab yab ' ^ab — 0 

“ to yah ■rab>0 


( 11 ) 
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Regarding the energy evolution an additional signal velocity 




u is adopted, and the viscosity term is defined as 


du\ 

dtjdi.s,a 


= -E 


m 

Pab 


^0:V^{g{\ab ■ + ^u^sig,u^ab 


rab-'^aWab 

( 12 ) 


where the signal velo city for the energy ju mps is chosen to be 

^sig,u Wab 


fail as in 
In order to proper 


Wadsley et al. (2008 1 . 
y compare these two different approaches 


and derive a Shakura & Sunyaev-like viscosity value for the dif¬ 
ferent simulations, we adopted the relations given by |Meru &| 
[Bat^ ( |20T2l l 


I h 3 „//! 


for the [Monaghan & Gingold| ( [l983| l formalism, and 
31 /i 9 


(13) 


(14) 


for the |Monaghan| ( [TW7| l one, where h is the averaged smoothing 
length, and H is the disk scale height. 

The initial average value of ass for our vine run is 0.02, sim¬ 
ilar to that in [Xiang-Gruess & Papaloizou] ( |2014| ). About twice 
larger is the ass in the phantom run with an average initial value 
of 0.04. 


Initial conditions — The scenario we have modelled includes a 
binary system composed of 2 equal mass stars iWp - m^- \ Mq. 
One of them (which we will call primary) harbors a protoplan¬ 
etary disk with a Jupiter mass planet Mp = 1 Mj embedded in 
it on an orbit with initial semimajor axis cip = 5 au. Different 
inclinations between the two components of the binary has been 
studied, in a reference frame centered onto the primary star. The 
disk is modelled with 600000 SPH particles, about 3 times the 
number used by Xiang-Gruess & Papaloizou ( 2014| l, and it ex¬ 
tends from 0.5 to 30 au, with a surface density profile 


E = Eor 


- 0.5 


(15) 


where Eo is defined such as the total disk mass is Md = 0.01 Mq. 
The choice of the initial parameters are such that the disk 


will precess as a solid body. In fact, according to Papaloizou & 
|Terquem] ( |1995| ) and jLarwood et al.] ( |1996| l, the condition for this 
behavior is 


Pout/Pd < H/R 


(16) 


In our initial setup this condition is fully satisfied since H/R 
is 0.037 while the ratio between Wp/Qa is approximately 0.01. 


This is confirmed also by the calculations of Xiang-Gruess & 
[Papaloizou] ( |2013| l. As reference frame we adopt the initial plane 
of the disk and planet orbit while the binary orbit has different 
initial inclinations defined respect to this plane. The initial semi¬ 
major axis of the companion star is ab = 100 au and it is set on a 
circular orbit. 


4. Initial mutual inclination of 45° 

The model where the inclination between the binary orbit and the 
initial disk-nplanet plane is 45° is assumed as standard model. 
The simulation with vine was halted after 11 600 orbits of the 
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planet (125 000 yr) while that with phantom was stopped after 
about half of that period (60000yr). The run with vine required 
approximately 8 months of CPU on a 24 processor machine 
while that with phantom 4 months on a 32 processor machine. In 
both simulations we find that in the first 6 000 yr the behavior is 
similar to that described in Xiang-Gruess & Papaloizou] ( [2014 1 
but, soon after, the dynamical evolution substantially changes. 
The inclination of the planet ip, that initially follows closely that 
of the disk id, decisively departs from id and grows at a lower 
rate. In Fig. |^both ip and id are shown as a function of time 
for the run with vine and phantom. While the periodic oscilla¬ 
tions of id continue with an almost unaltered frequency of about 
2 X 10"^ yr, the planet inclination ip grows on a much longer 
timescale reaching a maximum after about 8 x 10“* yr. There are 
clear indications of interaction between the planet and the disk 
in particular when the two inclinations are comparable like after 
2.5 X 10^, 5 X 10^ etc. At these times the planet crosses the disk 
and it interacts with its gas particles. 

There are some minor differences between the vine and the 
PHANTOM models possibly related to the different handling of vis¬ 
cosity and accretion onto the planet. A circumplanetary disk de¬ 
velops in the initial phase of the run with vine when ip ~ id 
but later on when the planet departs from the disk and periodi¬ 
cally crosses it with high mutual inclination it is dissipated. In 
both simulations dynamical behavior of the planet shows a clear 
competition between the gravitational force of the binary and the 
interaction with the disk, but the binary secular perturbations fi¬ 
nally dominate. The inclination of the planet roughly follows a 
pure N-body trend, shown as a dotted light-blue line in Fig. ]^ 
even if ip is strongly perturbed by the disk. The detachment be¬ 
tween the planet and the disk can be seen also in Fig. [Tj where a 
3D rendering of the system is displayed at an optical depth of ~ 
2 . 


In Fig. J^both simulations show a slow damping of the incli¬ 
nation oscillations, a behavior already suggested by[Martin et al.[ 
( ]20T4) . 

As discussed in Sec.]^ the distinct evolution of the planet and 
disk is related to their different precession timescales around the 
binary axis which translates in the inclination evolution shown 
in Fig.]^ As an additional indication of the different evolution of 
the planet respect to the disk in Fig. ]^ we illustrate the evolution 
of the precession angle of the disk Ld and planet Lp angular 
momentum around that of the binary star Lb, calculated respect 
to the initial orbital plane ( jXiang-Gruess & Papaloizou|[2014] 
[Larwood et al.|1996] l 


COSySp = 


Ld X Lb 

|Ld X Lb I 


(17) 


where u is the fixed unit reference vector for the initial orbital 
plane. The evolution timescale of the planet precession rate is ~ 
10 times slower than that of the disk, confirming the theoretical 
prediction. 

To test whether a potential warping of the disk due to binary 
perturbations is fully responsible for the decoupling of the planet 
from the disk in our models, as described in ( |Terquem|2013[l, we 
checked the disk shape every 300 planetary orbits. In Fig. [^ we 
show the disk evolution at different times and there is no evi¬ 
dence of a marked warping. In addition, we also performed a 
simulation where a smaller disk was considered. In Fig. ]^ the 
PHANTOM code was used to model the evolution of a disk with 
the same mass as our standard disk but extending only to 15 au in 
radial distance. The warping for such a disk, being less affected 
by the binary perturbations extending well within the tidal trun¬ 
cation radius ([Artymowicz & Lubow[[T994]l, is expected to be 
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Fig. 1. 3D rendering of the disk + planet system in the vine run. The 
plot is created through a ray tracing process, from an observer point of 
view placed at ~ 150 au from the disk. The disk is shown at an optical 
depth of ~ 2, defined as the cross section per SPFl particle mass. 



Time ( lO'^ yr) 


Fig. 2. Evolution of the disk and planet inclination respect to the initial 
disk+planet plane in the two runs with vine and phantom, respectively. 
The dashed magenta (vwe) and green (phantom) lines show the planet 
inclination ip while the continuous black (vine) and red (phantom) line 
illustrates the disk inclination ij. After about 6000yr, the planet evo¬ 
lution departs from that of the disk in both simulations and ip grows at 
a slower rate compared to ij. The light blue dashed-dotted line illus¬ 
trates the evolution of the planet inclination in a pure N-Body problem 
Star-Planet-Companion star. 



negligible. Even in this case the planet decouple from the disk 
in a way very similar to that shown in Fig. even if on slightly 
different timescales due to the different disk configuration and 
density. Initially, the smaller disk has a superficial density that is 
about 6 times higher compared to our standard disk. As a con¬ 
sequence, when the planet decouples from the disk, the repeated 
crossing of the disk plane lead to a different dynamical evolution 
since the frictional force depends on the local disk density. The 
increase of only 25% in the precession period of the disk initially 
extending over 15 au, instead of the expected factor 2 on the ba¬ 
sis of eq. 3), is due to the radial spreading of the disk due to the 
strong binary perturbations. The standard disk with initial outer 
radius of about 30 au cannot spread beyond the tidal truncation 
limit and, as a consequence, its outer radius remains constant. 
On the other hand, the smaller disk has room to expand in re¬ 
sponse to the binary gravitational perturbations extending after 
one cycle (about 2 x 10^* yr) beyond 25 au. The outcome of the 


Fig. 3. Evolution of the disk and planet precession angle respect to 
the initial disk-l-planet plane, as defined in eq. \ni The dashed green 
line shows the planet precession angle /ip while the continuous red line 
illustrates the disk’s one/ij. 


numerical simulations rule out warping as responsible for the de¬ 
coupling of a planet from its disk in our dynamical configuration 
since, according to ( |Terquem|2013| l, a severe warping is needed 
to affect the planet inclination. However, this mechanism may be 
an important one for closer binaries, increasing the decoupling 
rate between planet and disk. 

The planet migration, in Fig. [6] is significantly slower com¬ 
pared to that found by (see [Xiang-Gruess & Papaloizou |2014| 
Fig. 3). Apart from the initial fast migration rate when the planet 
is still embedded in the disk, when the planet detaches from the 
disk plane it is the friction developing during the periodic cross¬ 
ing of the disk by the planet that dominates the semimajor axis 
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Fig. 4. Disk shape shown every 300 planetary orbits when the inclina¬ 
tion of the planet is decoupling from that of the planet. The disk warping 
is negligible. 



Fig. 5. Comparison of the planet and disk inclination evolution in the 
standard case, where the disk extends to 30 au, and a model with a 
smaller disk, extending to 15 au. No warping is expected in the latter 
case but the planet decouples from the disk as in our standard case with 
a larger disk. 


evolution and migration. According to Teyssandier et al. ( 2013| l, 
the planet experiences a significant friction while crossing the 
disk due to aerodynamic drag and to a dynamical drag due to the 
gravitational scattering of the disk particles. The latter dominates 
since the ratio between the two friction forces can be expressed 
as ([Teyssandier et al.|2013|l 


^ aer 
F dyn 




( 18 ) 


where Mp and Rp are the mass and radius of the planet, respec¬ 
tively, and a is its semimajor axis. In our initial conhguration 
where a = 5 au, H{r) - h^a ~ 0.19 and assuming a radius of 
the planet Rp ~ IRj ~ 0.001 au we get a ratio lower than 0.02. 
In this situation the drag force is mostly due to the particle scat¬ 
tering during the planet crossing. There may also be an angular 
momentum exchange between the planet and the disk when the 
planet approaches the disk, but it has probably less effect on the 
long term drift. 

The migration rate appears to be more irregular in the run 
with VINE respect to that with phantom. This different behavior is 
possibly due to the mass accretion onto the planet in the phan¬ 
tom run (the planet is modelled as a sink particle) which it is 
not implemented in the vine run. In this last the planet cannot 
growth and gas is temporarily trapped and removed around it 
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during the disk crossing possibly leading to a noisy migration. 
In the run with phantom, the planet has reached a mass of 1 .6 Mj 
after about 4 x 10"^ yr and then the mass growth is reduced. 



Fig. 6. Evolution of the semimajor axis of the planet in the vine and 
PHANTOM runs, respectively. The planet migrates because of the friction 
with the gas when it periodically crosses the disk plane. 


5. Initial mutual inclination of 60° 


Even in the scenario with is - 60° we have run two different 
simulations, one with vine and one with phantom. In Fig. |7]we 
show the evolution of the disk and planet inclinations over a 
shorter timescale compared to the previous simulations where 
/b = 45°. However, it is already clear that the planet abandons 
the plane of the disk and evolves independently under the action 
of the companion star in both simulations. The disk influence 
appears to be a perturbing force also in this different scenario. 
The planet inclination, as in the /b = 45°, follows that of the 
disk only for a few 10^ yr after which the planet detaches from 
the disk following an N-body behavior while the disk perturbs 
the dynamics especially when the planet crosses the disk. In the 
two models the planet evolution is different and this can possi¬ 
bly be ascribed to the diverse way of handling the gas trapped 
around the planet and to the different behavior of the disks in 
the two cases. After an initial quick increase of the disk inclina¬ 
tion, the growth halts at about - 90° and damped oscillations 
are thereinafter observed in both cases. However, the increase 
in inclination appears more irregular in the run with vine and 
stronger damped oscillations are observed immediately after the 
hrst maximum. In the run with phantom, the evolution appears 
more regular and the damping more progressive. This different 
behavior may be ascribed to the different viscosity which seems 
to play a more relevant role on the evolution of the disk inclina¬ 
tion when i'b is higher. 

It is noteworthy that the damped oscillations recall, at some 
extent, those observed in jMartin et al.| ( ?014| l. However, the pa¬ 
rameters of our model significantly differ from theirs so the ex¬ 
act same behavior is not expected. In particular, [Martin et al.j 
( 2014[ l use 1) different temperature and density profiles, 2) their 
viscosity parameters are set in order to have a constant value of 
ass all over the disk, and 3) their disk is ten times less massive 
compared to ours. As a consequence, the evolution of the disk 
inclination is only qualitatively similar, but it is interesting that 
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in both scenarios damped oscillations are observed. An initial 
hint of damped oscillations of the disk inclination is also seen in 
the model with - 45° (Fig.[^ even if, as already observed, it 
appears feeble. 

A potential contribution to the strong damping in the - 
60° case is possibly related to the exchange of disk mass between 
the two stars. After about 2 x 10"^ yr from the beginning of the 
simulation the secondary star has gained a gaseous disk whose 
mass is about 5% that of the primary disk. The mass exchange 
is significantly lower in the - 45° case where the acquired 
disk has a mass of only 0.05% of the initial disk mass. As a 
consequence, the viscosity of the disk plays a role also in an 
indirect way by forcing mass exchange for higher values of 
leading to a damping of the disk inclination oscillations. 



Time ( lO^yr) 


Fig. 7. Evolution of the disk and planet inclination respect to the initial 
disk+planet plane in the case with Ib = 60° in the vine and phantom 
runs, respectively. The symbols are the same as in Fig.|^ 



Fig. 8. Evolution with time of the disk in the phantom model. Some 
mass leaves the primary disk and forms a ring around the secondary star. 
This behaviour was already observed in |Picogna & Marzari| ( |2013t . No 
significant warping is observed. 


Even in this case we tested the potential warping of the disk 
by plotting the disk evolution every 300 orbits of the planet. In 
this more inclined configuration we observe some mass transfer 
from the disk around the primary to the secondary star where 
a ring of gas slowly builds up. It is interesting to note that the 
disk eccentricity is higher compared to the case where is = 45" 
and comparable to that observed in Martin et al. ( |2014 1 . This is 
why the mass flux from the primary towards the secondary star is 
more marked. In any case, the amount of mass extending out of 
the primary disk is not enough to cause a significant gravitational 
perturbation on the planet. 


6. Discussion and Conciusions 


We have shown that a Jupiter-like planet embedded in a circum- 
stellar disk will evolve almost independently from the disk in 
response to the perturbations of a misaligned binary stellar com¬ 
panion for values of the semimajor axis as around 100 au. For in¬ 
clinations of the binary plane of the order of 45° and 60°, after an 
initial coupled evolution, the planet dynamically detaches from 
the disk and its inclination oscillates on a significantly longer 
timescale compared to that of the disk. This occurs since the N- 
body secular perturbations of the companion star dominate over 
the damping force of the disk which tends to drag the planet back 
into its median plane. Our simulations have 3-times the resolu¬ 
tion of those performed by Xiang-Gruess & Papaloizou ( 2014) l 
and were performed with 2 different up to date SPH codes, vine 

and PHANTOM. 

Our results show that one of the fundamental requirements 
of the model of |B aty gm] ( |20 1 2) l to explain the observed spin- 
orbit misalignment of some exoplanets is not met for small bi¬ 
nary separations and a reasonable choice of physical parameters 
of the disk and planet. However, the migration of the planet ap¬ 
pears to occur at a significant rate even if the orbit of the planet 
and the disk plane are not aligned. This is due to the friction with 
the gas during the repeated crossing of the disk by the planet. As 
a consequence, the inward drift of the planet would not be due 
to type I/II migration, as argued for wider separations by |Crida| 
|& Batygm ( |2014[ ), but to the dynamical friction with the disk. 
This would possibly lead to conclusions similar to those given in 
|Batygin| ( |2012| l concerning the long term evolution of the planet 
spin axis and semimajor axis migration, but the physical mecha¬ 
nism would be different. The formation of hot Jupiters with spin 
misaligned respect to the stellar equator would in this case not 
occur because the planet follows the disk in its precession, but 
by a combination of N-body Kozai cycles and repeated crossing 
of the disk plane by the inclined planet evolving out of the disk. 
Even if the disk oscillations are damped over a long timescale to 
a given value of inclination, the continuous crossing of the disk 
by the planet following the Kozai-Lidov cycles would grant a 
significant migration. 

Of course, a deeper exploration of the parameter space is 
needed to have a more complete scenario of the behavior of the 
system. We have also explored only one particular step of the 
system evolution while previous steps have also to be investi¬ 
gated. Ad example, when does the planet begin to evolve inde¬ 
pendently from its birth disk? At which stage of its growth? For 
which value of mass? 

We have shown in our simulations that a Jupiter mass planet 
exits almost immediately from the disk, but does this occur also 
for smaller mass planets? Possibly not, since the final infall of 
gas on the planet core is a fast process occurring on a timescale 
shorter than the decoupling between the planet and disk inclina¬ 
tion. However, the problem of growing a planet core in a disk 
perturbed by an inclined binary is a critical process. The con¬ 
densation of dust into larger pebbles and, eventually, kilometer¬ 
sized planetesimals may be strongly perturbed in particular when 
the planetesimals decouple from the gas. For small inclinations 
between the circumprimary gas disk and the binary orbital plane 
|Xie & Zhou| ( |200^ and |Xie et al.| ( |201^ showed that planetesi- 
mal accretion may be fast because size segregation among plan¬ 
etesimals favors low velocity impacts. For larger inclinations, 
the nodal longitude randomization forced by the companion star 
may lead to the dispersion of the planetesimal disk reducing the 


chances of planet formation. However, Marzari et al. ( 2009bl l 
showed that for binary semimajor axes larger than 70 au the 
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nodal longitude randomization timescale is longer than the plan- 
etesimal accretion process and core formation might occur. As a 
consequence, planet formation might occur within the disk even 
in presence of an inclined binary until the planet reaches a Jupiter 
mass or even more. After that, it would decouple from the disk 
and follow and independent evolution. 

Additional simulations are needed to better explore the prob¬ 
lem and different initial values have to be sampled for the planet 
mass, the binary orbital elements, and the disk parameters (like 
mass, viscosity, temperature distribution, density profile and so 
on) which can influence the propagation of perturbations. Run¬ 
ning more numerical models may help to understand not only the 
response on the long term of the disk to the binary perturbations 
and the possible occurrence of damping, but also the dynamical 
evolution of the planet. The problem with these type of simula¬ 
tions is the large amount of CPU required by any run of an SPH 
code in this configuration. 

Even if the migration of the planet in this scenario is due to 
the friction with the the gas of the disk, it remains an open prob¬ 
lem how to drive the planet very close to the star. Isothermal 
simulations of disks in binaries (Kley et al.|2008[|Marzari et al. 


|2009a| l show on the long term the formation of an elliptical hole 
in the density distribution close to the primary star. This would 
evidently halt the planet migration once the planet moves within 
the hole where either type I/II migration or, as we suggest in this 
paper, friction with the gas during the crossing of the disk would 
be switched off. However, more recent simulations with radiative 
disks do not show the formation of such an internal low density 
region ([Mliller & Kley 2012 [Marzari et al. 2012[ Picogna & 
|Marzari|2013| ). This is an additional aspect that should be inves¬ 
tigated to really understand the dynamical evolution of inclined 
hot Jupiters in misaligned binary star systems. 
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